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ABSTRACT 


Wind fields were measured with the ground-based NASA Marshall Space 
Flight Center (NASA/MSFC), Huntsville, Alabama, lidar, NOAA Wave Propaga- 
tio l Laboratory (NOAA/WPL), Boulder, Colorado, and with the NASA B-57B 
instrumented aircraft. The remotely sensed winds are compared with the 
in situ aircraft measurements. Three flight plans were carried out 
during the two different field programs. At NASA/MSFC the aircraft 
circled while the lidar scanned conically on May 10, 1983, and the air- 
craft flew 6° approach path al:ng the fixed lidar beam on May 12, 1983. 

In Boulder, Colorado, on February 7 and 9, 1984, the aircraft flew an 
approach along the lidar beam directed south-north (parallel to the 
mountain range) and a climbout along the lidar beam which alternately 
shifted east-west (perpendicular to the mountain range). Turbulence 
intensities and spectra were calculated from the temporal fluctuations 
in the lidar-measured radial wind speed component. It should be noted 
that time histories of the lidar wind represent values spatially averaged 
over a 300 m volume element. The lidar winds were sampled at approxi- 
mately 2 times per second whereas the aircraft measurements were sampled 
at 40 times per second. 

The second moment or Doppler frequency spectra width of the lidar 
measurements was also compared with turbulence intensities measured by 
the aircraft. Variable sample sizes of pulses were averaged in order to 
resolve turbulence information from the lidar spectra width. 

It is concluded that these field tests provided unique sets of data 
to examine the mean wind and turbulence measurements made by remote 
sensing instruments. The comparison of aircraft-measured turbulence 
intensities and spectra with lidar time histories of radial wind speed 
were in good agreement. Although the magnitude of lidar second moment 
(or spectral width) is 4 to 5 times higher than turbulence intensity, 
variation cf the second moment does contain information representative 
of the actual measured turbulence. 
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1.0 INTRODUCTION 


A technique of remotely sensed wind measurement was described by 
Huffaker et al. (1970) by applying the concept that velocity can be 
determined frrm the Doppler shift in light scattered by particles in the 
atmosphere (Bi'bro 1980). A comparative study between a cup anemometer 
and a continuous wave C0£ Laser Doppler Velocimetry (LDV) has been per- 
formed by Lawrence et al. (1972). Also, Brashears and Hallock (19/6) 
reported the first wind profile measurement using a pulsed system. The 
measurements were compared with National Weather Service (NWS) wind 
soundings. The term "lidar" is a general term of this new remotely 
sensing system and is an acronym derived from light detection and 
ranging. One substantial improvement of wind measurement using pulsed 
lidar was made with the appearance of a transverse-excited atmcshperic 
(TEA) laser. Hall et al . (1984) reported the comparison of wind mea- 
surements using a TEA configuration lidar system and tower anemometer, 
rawinsonde, and the profiler. However, most comparison studies utilized 
a conical scanning technique known as velocity azimuth display (VAD) to 
resolve three-dimensional winds by using a single l^dar system (Cliff 
and Huffaker 1974). Since the velocities determined by the VAD technique 
are the result of large-scale spatial averaging, the information on • 
small-scale turbulence is lost. The lidar-measured wind and turbulence 
at different altitudes (range gates) have not been examined before by an 
in situ measurement, such as available from an instrumented aircraft. 

Two field tests with different comparison configurations were carried 
out in this regard. 

The first field test was conducted at NASA Marshall Space Flight 
Center (NASA/MSFC) on May 10 and 12, 1983, in order to compare NASA/MSFC 
lidar measurements with NASA B-57B aircraft measurements. On May 10, 
lidar-measured winds were determined using a conical scan while the 
aircraft flew circular flight paths at several altitudes. On May 12, 
the lidar was fixed at a 6° elevation angle while the aircraft flew 
approach paths roughly parallel to the lidar beam. Results of the May 
12 test have been presented in a previous report (Frost and Huang 1983). 


1 



The second field tests were conducted at Boulder, Colorado, on February 
7 and 9, 1984. Similar to the May 12 tests, the aircraft (NASA B-57B) 
flew an approach along the NOAA/WPL lidar beam directed south-north 
(parallel to the mountain range) and a climbout along the Ildar beam 
which shifted to an east-west (perpendicular to the mountain range) 
flight path. 

T he mean wind fields, the turbulence intensities, and the turbu- 
lence spectra determined from measurements by both systems are in good 
agreement. Turbulence intensities and spectra were calculated from the 
fluctuations with time for the 1 iuar-measured radial wind speed compo- 
nent. The second moment or Doppler frequency spectra width of the lidar 
measurements was also compared with turbulence intensities measured by 
the aircraft. Variable sample sizes of pulses were averaged in order to 
resolve turbulence information from the lidar spectra width. 

In this report, the instrumentation involved in obtaining the data 
and an outline of the field test plan is given in Section 2.0. In 
Section 3.0, the method of comparison for each of the field tests is 
described and the comparison and interpretation of the results from both 
measurements ( i . e . , lidar and aircraft) are also presented. 

It is concluded that these field tests provided unique sets of data 
to examine the mean wind and turbulence measurements made by remote 
sensing instruments. The comparison of aircraft-measured turbulence 
intensities and spectra with lidar time histories of radial wind speed 
was good agreement. It is also concluded that the lidar second moment 
or Doppler lidar spectra width holds promise for being a turbulence 
indicator. 
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2.0 EXPERIMENT 


2 . 1 Instrumentation 

As mentioned previously, two field tests were carried out for this 
study using different Ildar systems for each test. During the first 
field test the instrumentation consisted of the NASA/MSFC C0 2 Doppler 
Ildar and the NASA B-57B instrumented aircraft and was conducted at 
NASA/MSFC, Huntsville, Alabama. The second field test was carried out 
at Boulder, Colorado, where the NOAA'WPL lidar system was operated to 
compare with NASA B-57B measurements. The Instrumentation used In the 
first field test is discussed in a previous report (Frost and Huang 
1983). Details of the NASA/MSFC Doppler lidar are also given by Bilbro 
and Vaughan (1978), Jeffreys and Bilbro (1975), and Lee (1982). Infor- 
mation on NASA's B-57B aircraft is described by Camp et al. (1983), 
Campbell et al. (1983), and Theon (1985). Similar to the NASA/MSFC 
lida>, the NOAA/WPL lidar system is also a C0 2 pulse Doppler li The 
difference in operation between the two lidar systems are listed in 
Table 2.1. Details of the NSAA/WPL lidar system is provided by Post et 
al. (1981). The main difference between the systems is their configura- 
tion. The configuration of the NASA/MSFC pulse lidar is a master oscil- 
lator power amplifier (MOPA), while NOAA/WPL pulse lidar is a hybrid, 
transverse-excited atmosphere (TEA) configuration. One shortcoming of 
the MOPA configuration is the relatively low per-pulse energy, typically 
10 to 30 mj (Bilbro 1980). It takes at least 50 pulses to obtain a 
meaningful velocity measurement. Moreover, it is poor as a second 
moment (spectral) estimator, due to the signal -to-noise ratio (SNR) 
being low. These deficiencies have been improved substantially by the 
TEA laser system. The TEA laser which is used at NOAA/WPL, produces 100 
mJ of energy at a 12 Hz rate, which provides sufficient frequency 
stability to allow wind velocity measurement, and also to give meaning- 
ful output as a second moment estimator. 

2.2 Field Test Design 

This report emphasizes the data from the NASA/MSFC May 10, 1983, 
test and the NOAA/WPL field tests of February 7 and 9, 1984. During the 
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TABLE 2.1. Comparison of the Lidar System Parameters Between NASA and 
NOAA Lidar. 


Parameter 

NASA/MSFC 

NOAA/WPL 

Source 

C\J 

o 

o 

co 2 

Wavelength 

10.6 ym 

10.6 ym 

Configuration 

MOPA 

TEA 

Pulse Energy 

25 mJ 

100 mJ 

Pulse Duration 

2 ys 

2 ys 

Pulse Per Second 

110 pps 

12 pps 

Number of Pulse Averaged 

50 

6, 24, or 48 
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field test of May 10, 1983, the Doppler Ildar was operated In a conical 
scanning mode, or the velocity azimuth display (VAD) mode. Scans were 
carried out at elevation angles of 6°, 9°, 13°, 19°, 26°, and 32°. The 
aircraft flew circular flight paths at Increzslng altitudes in order to 
approximately capture the Ildar beam as Illustrated in Figure 2.1. 

Figure 2.2 shews one circle of a VAD at the 6° elevation angle. Each 
vector shown in Figure 2.2 represents the radial component of wind 
velocity. The mean wind direction can be detected as roughly 200 ±20° 
from the plot. Wind velocity data of four different runs (Runs 13, 15, 
18, and 20) from aircraft measurements are also available. Table 2.2 
lists the altitudes and mean radii for the circular patterns nf the 
aircraft trajectories. Based on the aircraft altitude, the range gate 
for the corresponding altitude was selected ior each run in order to 
compare the two measurements. Figure 2.3 shows the horizontal wind 
vector of the aircraft measurement for Run 13 at an altitude of 592 m 
above mean sea level (MSI). The mean wind direction Is also in the 
range of 200° to 240°. The NASA/MSFC tower measurements (Figure 2.4), 
which were obtained just before taking the measurements used in the 
lidar and aircraft comparison study, also shows similar direction at 187 
m MSL (20 m above local terrain). The figure clearly Indicated that the 
mean wind direction approaches 200° near the end of the time period over 
which the measurement was made. This is in agreement with both lidar 
and aircraft measurements at higher altitudes. 

The second field test was conducted at Boulder, Colorado. The 
N0AA/WPL lidar was located on the northeastern corner of Table Mountain. 
On both test days (i.e., February 7 and 9, 1984), the lidar was fixed at 
a 4.5° elevation angle for both the 20C° and 290° azimuthal from true 
north. The aircraft flew an approach (slope aporoxima y 4.5°) into 
the lidar beam directed at 200° azimuth (roughly parallel to the moun- 
tain range) and a climbout along the lidar beam which shifted to an 
azimuth of 290° (roughly perpendicular to the mountain range). Figure 
2.5 shows the lidar beam location relative to the terrain. The synoptic 
flow patterns of both test days are different. The prevailing wind 
direction was 30° from true north (parallel to the mountain range) on 


5 





Distance (m) 


OBKtiNAL PAG- .S 
bfifOOR QUALITY 



Distance (m) 

Figure 2.2. Radial wind vectors of lidar measurement from range gates 6 
to 30 scanning 360° at 6° elevation angle. 
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TABLE 2.2. Comparison of B-57B and Lidar Data. 


Flight 40 
Run No. 

B-57B Mean 
Radius 
(m) 

B-57B 
Altitude 
MSL (m) 

Time 

(sec) 

Lidar 
Elevation 
Anqle (deq) 

Lidar 
Range 
Gate No. 

_idar Mean 
Radius (m) 

13 

4,400 

592 

350 

6 

14 

3,920 

1 5 

4,500 

913 

371 

9 

17 

4,628 

18 

5,120 

1,226 

340 

13 

16 

4,530 

20 

5,100 

1 ,834 

353 

19 

17 

4,803 
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Figure 2.3. Horizontal winds along the B-57B flight pattern; Run 
May 10, 1983, at NASA/' MS FC . 
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Figure 2.5. Two fixed lidar beams (200° and 290° azimuth) relative to the 
contour at Boulder, Colorado. 
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February 7 while it was 310° (oblique relative to the mountain range) on 
February 9. For the flow oblique relative to the mountain range, the 
general flow pattern is more like a vortex shed from a ridge of the 
mountain range. This is similar to a flow approaching a fence at an 
oblique angle. Figure 2.6 schematically illustrates a vortex flow 
pattern while Figure 2.7 shows how laboratory wind tunnel simulations 
have identified vortex flow off the Rock of Gibraltar (Cook et al. 

1978). 

Figure 2.8 shows the time histories of radial wind component which 
was measured at 290° azimuth on February 9. It clearly shows the reverse 
flow associated with the vortex flow at lower altitudes. Detailed 
analysis of the data and comparison of both measurements will be dis- 
cussed in Section 3.0. 
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Figure 2.6. Illustration of leading edge vortex due to oblique 
flow over a long fence or ridge. 



Figure 2.7. Wind tunnel observation of the vortex for wind from 200° 
true north at the Rock of Gibraltar (Cook et al. 1978). 
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.8. Time history of 1 idar-measured wind vector at 290 c azimuth, 
4.5° elevation relative to the terrain on February 9, 1984, 
Boulder, Colorado. 
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3.0 COMPARISON OF LIDAR MEASUREMENTS WITH 
AIRCRAFT MEASUREMENTS 


Since the study was to compare Doppler lidar-measured winds and 
turbulence with in situ aircraft measurements and since the instrument 
orientations of both measurements are different, the aircraft data were 
reconstructed according to the Ildar beam configuration for each field 
test. It is also necessary to assure that in situ measurements were 
compared with lidar measurements at corresponding time and spatial 
positions. In this section, results from the May 10 data will be 
discussed first. An investigation of correlation between lidar ampli- 
tude and turbulence intensities is then given, and finally results for 
the February 1984 field test are presented. 

3.1 NASA/MSFC Field Test 

As indicated in Table 2.2, it takes 6 minutes for the aircraft to 
complete a 9000 m diameter circle. Figure 3.1 shows the one-second 
average horizontal wind vectors measured by the aircraft. A total of 
four runs at different altitudes are presented in this figure. The 
figure shows the spatial wind variation at each level and also a mean 
wind velocity increase with altitude. These wind vectors were plotted 
according to the longitude and latitude coordinates which were recorded 
on the flight data tape. An error in the aircraft's inertial navigation 
system (INS) was encountered during the flight tests. The aircraft's 
longitude and latitude data indicate the lidar site is not enclosed in 
the circular pattern, which in fact is not the case. Evidence of INS 
problems was also shown for the May 12 test (Frost and Huang 1983). In 
order to compare with the lidar measurements, these circular patterns 
were shifted such that the centers of these circles were at the lidar 
site and in accord with visual observations. As shown in Figure 2.2, 
the lidar-measured wind is the radial component along each range gate 
for given elevation and azimuth angles. However, the recorded aircraft 
measurements are the three components (east-west wind, W £ , north-south 
wind, W^, and vertical wind, W^) of wind velocity at the aircraft's 
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position In the earth coordinate system, as shown in Figure 3.2. The 
radial component of the aircraft measurement, W R , was computed according 
to the following equation: 

0 

W R = (W E sin <t> + cos 41) cos 0 + sin e 

where <}> is the azimuth angle and 0 Is the elevation angle. Both <t> and 0 
can be calculated using the longitude and latitude data of the aircraft 
relative to the coordinates of the lidar site. Radial components of 
lidar measured wind at the range gate which was located at the aircraft 
altitude were then selected for comparison with the in situ measurements. 
The corresponding gate numbers are given in Table 2.2. 

Figure 3.3 shows a comparison of mean radial wind speed from the 
Ildar and aircraft measurements. Since the sampling rate of the aircraft 
measurements is 40 Hz, and it takes 6 minutes to complete one circle, the 
output data for one run is roughly 17,000 data points. The sample rate 
of the lidar is 2 Hz, and it takes 75 seconds to complete one circle, 
thus 150 data points were collected for each 360° scan. Although both 
measurements are very nearly at the same altitude, the sampling volumes 
could be separated by 100 m to 500 m as estimated for the circle radii 
(see Table 2.2). These effects may explain the difference between the 
two measurements. Generally speaking, however, agreement between the 
aircraft and the lidar radial mean velocity components is good. Both 
measurements show the same range of mean wind direction between 200° to 
240°. The spike seen in the lidar measurement in these figures is believed 
to be caused by a reflection of the lidar beam from a meteorological 
tower near the lidar site. 

Since the lidar scanned conically, the radial component of wind 
velocity appears to have a sinusoidal shape as shown in Figure 3.3. For 
a given altitude (or range gate), if the wind direction and speed a*"° 
constant, the radial wind component must be a perfect sinusoidal curve. 

A sine curve was fit to each data set. The difference between the 
radial mean wind and the sinusoidal curve fit is taken as the radial 
turbulence component. Figure 3.4 shows the original radial wind velocity 
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the aircraft-measured wind 



Wind Speed (m/s) Wind Speed (m/s) 




Azimuth Angle 
(b) Run 15, h = 913 m MSL 

Figure 3.3. Comparison of radial wind speeds between aircraft and lidar 
measurements . 
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Wind Speed (m/s) Wind Speed (m/s) 
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Radial Turbulence Radial Mean Wind Radial Turbulence Radial Mean Wind 

(m/s) (m/s) 


Time (sec) 

(a) h = 590 m MSL , 6 = 6° 



0 . 00 . 160 . 210 . 320 . 100 . 100 . 560 . 610 . 720 . 000 . 


Time (sec) 

(b) h = 942 m MSL, e = 9° 


Figure 3.4. Time histories of 1 idar-measured wind velocity and its 
turbulence component at different runs. 
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0 . 00 . 160 . 210 . 320 . 100 . 100 . 560 . 610 . 720 . 

Time (sec) 

(d) h = 1377 m MSL, e = 19° 
Figure 3.4. (continued). 
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and the turbulence component after a sinusoidal curve fit has been 
removed. 

Similarly the radial turbulence component of the aircraft measure- 
ments can be extracted from the radial wind velocity as shown in Figure 
3.5. Since the aircraft flew only one circle for each run, only one 
period of a sinusoidal curve needs to be removed. With the time histories 
of the turbulence component resolved in the preceding manner, the turbu- 
lence intensity can be calculated. Figure 3.6 shows the turbulence inten- 
sity at different altitudes along with the mean radial velocity for given 
mean wind directions for both measurements. The wind speed profile 
contains 2 to 3 m/s difference, while the turbulence intensity shows 0.3 
to 0.5 m/s difference between the aircraft and lidar measurements. The 
fact that the magnitude of aircraft-measured wind is 2 to 3 m/s higher 
than lidar-measured wind most of the time might be partially attributed 
to the Schuler frequency of the aircraft navigation system. A detailed 
investigaton of this is recommended. The comparison of lidar second 
moment (spectral width) with aircraft-measured turbulence intensity is 
not possible since the lidar second moment data is not meaningful at 
the altitude for which aircraft measurements were made. 

In order to investigate the turbulence measurements further, the 
turbulence energy spectra were computed. The spectra from aircraft- 
measured turbulence were computed by averaging ten segments of the time 
record, each segment contained 1024 data points. The spectra of lidar- 
measured turbulence, on the other hand, were computed by averaging 4 t( 8 
segments of the time history, each segment containing a total of 128 
points. The results are shown in Figure 3.7. Since different sampling 
rates are involved, the spectra distributions of lidar-measured turbulence 
cover a range of 0.01 Hz to 1 Hz, while spectra of aircraft-measured 
turbulence cover from 0.04 Hz to 20 Hz. The spectra distributions of 
both measurements agree well in the frequency band where they overlap, 
except at the highest measurement level, e = 19°. The deviation of the 
lidar turbulence spectra at this level from the aircraft turbulence 
spectra may be due to the length of lidar data time recorder being too 
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Radial Turbulence Radial Wind Velocity Radial Turbulence Radial Wind Velocity 
( m / s ) (m/s) ( m /s) (m/s) 



0. 10. 00. 120. 160 . 200. 210 . 200. 320 . 360. 100. 

Time (sec) 

(a) h = 592 m MSL, Run 13 



• io. l * ■ * * * * * 

0. 10. 00. 120. tbO. 200. 210. 200. 320. 360. 100. 

Time (sec) 

(b) h = 913 m MSL, Run 15 

Figure 3.5. Time histories of aircraft-measured radial wind velocity and 
its turbulence compo.ient at different runs. 
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Radial Turbulence Radial Wind Velocity Radial Turbulence Radial Wind Velocity 

(m/s) (m/s) (m/s) (m/s) 



0 . 10 . 00 . 120 . 160 . 200 . 210 . 200 . 320 . 360 . 100 

Time (sec) 

(c) h = 1226 m MSL, Run 18 



0 . 10 . 00 . 120 . 160 . 200 . 210 . 200 . 320 . 360 . 100 . 

Time (sec) 

(d) h = 1834 m MSL, Run 20 
Figure 3.5. (continued). 
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Figure j 6. Comparison of radial mean wind speed at mean wind direction 

and calculated turbulence intensity between lidar and aircraft 
measurements . 




JrtlGtNAL PiSK.:^ 

Oh POOR QUALITY 



f (Hz) f (Hz) 

(a) h = 592 m, e = 6° (b) h = 913 m, e = 9° 




0.01 0. IC 1.00 1C. 00 100.00 lCCC.CO 


f (Hz) 

(c) h = 1226 m, 0 = 13° 


f (Hz) 

(d) h * 1831 m, e = 19° 


Figure 3.7. Comparison of radial turbulence spectra computed from aircraft 
and lidar measurements. 


28 





short for this particular run (as shown in Figure 3.4d) or due to decreas- 
ing SNR at the higher range gates (#18). Figure 3.7 also shows that, 
generally speaking, the spectra distributions of both aircraft- and lidar- 
measured turbulence follow a -5/3 slope for the inertia subrange, which 
is a typical finding for atmospheric turbulence. 


3 . 2 Lidar Amplitude Rela te d to Turbulence Intensity 

As mentioned In the previous report (Frost and Huang 1983), the 
computed turbulence intensity for lidar-measured winds is extremely 
uncertain at the higher elevations (or higher range gates). Therefore, 
an effort was made in this study to determine if a calibration existed 
between turbulence intensity and lidar signal amplitude (or lidar signal 
intensity). Data from May 12, 1983, were utilized to investigate if such 
a relationship exists. Figure 3.8 is a plot of turbulence intensity 
versus lidar signal amplitude. 


The independent variable is the lidar amplitude as provided on the 
NASA/MSFC lidar data tapes. The dependent variable in Figure 3.8 is the 
ratio of turbulence intensity calculated from the time historip*; of the 
radial velocity from the lidar to that of the B-57B aircraft, respec- 
tively. Assuming that the B-57B measures the real turbulence intensity, 
a relationship between 1 idar-calculated turbulence and lidar amplitude 
would provide a correction factor for adjusting °|_-j cj ar * F ''9 ure 3.8 
suggests that a functional relationship of the form 0 Li ( jar^ 0 Aircraft = 
f(amp) may exist. Plotted in the figure are composites of eight runs. 

Data for each specific run is indicated by the number at the data posi- 
tions. The figure contains data from only range gates 9 through 18. 

Applying a curve fit routine, the calibration formula is found to be: 

f 0.11 (amp-50) 
e amp >_ 50 db 

f(amp) = 

1 amp < 50 db 

An inherent assumption, however, is that there is no variation in the 
aerosol content over the altitude interval investigated. If there is a 
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strong relation between Ildar amplitude and aerosol content. Figure 3.8 
may be simply Illustrating that phenomenon. In order to establish 
confidence In the potential relationship between Ildar amplitude and 
turbulence Intensity, further studies are required. 

3.3 NOAA/WPL Field Test 

The field test carried out in Boulder was similar to the May 12, 

1983, test at NASA/MSFC as Indicated In Section 2.2. The method used to 
analyze these data was the same as the method Indicated In a previous 
report (Frost and Huang 1983). The experiment was a two-day test, l.e., 
February 7 and 9, 1984. The aircraft flew a total of 16 paths. Including 
8 approaches along the 20° azimuth (parallel to 200° azimuth Ildar beam) 
and 8 climbouts along 290° azimuth for each day. However, only six 
runs of the February 7 flight and four runs of the February 9 flight 
overlapped with lidar measurement, as shown in Tables 3.1 and 3.2. 

Figures 3.9 and 3.10 show the wind vector for both days. It clearly 
shows that the prevailing wind on February 7 was blowing toward the 
mountain while it was blowing from the mountain on February 9. Some ol 
the flight trajectories appear to ue far from the lidar beam which again 
is believed to be the INS drift problem. The 1 idar-measured wind in 
each range gate is a spatial average over the 300 m length of the sample 
volume (2 ps pulse duration). The time history of the aircraft -measured 
wind is then averaged over a period corresponding to the length of time 
required for the aircraft to traverse the 300 m range gate along the 
flight path. As indicated in the previous report, two approaches to 
carrying out this averaging technique were adopted. One was to assume 
vertical homogeneity in the flow field as shown in Figure 3.11a and the 
other technique was to average the wind assuming horizontal homogeneity 
as shown in Figure 3.11b. 

Figures 3.12 and 3.13 compare the 1 idar-measured mean radial wind 
velocity with the aircraft-measured winds averaged over the corresponding 
300 m sample volume, assuming horizontal homogeneity at 20° azimuth arid 
290° azimuth, respectively. Figures 3.14 and 3.15 show the same comparison 
except assuming vertical homogeneity. It was difficult to find two 
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TABLE 3.1. Selected Runs of the February 7 Test. 



B-57B Aircraft Data 


NOAA Lidar Data 

Run 

Azimuth 

Sampling Time 
(MST) 

PRF 

Number of 
Pulse 

Sampling Time 
(MST) 

No. 

Angle 

Start to End 

(Hz) 

Averaqe 

Start to End 

2 

290 

11:46:42-11:49:19 

12 

6 

11:46:53-11:49:04 

3 

200 

11:56:42-12:00:27 

12 

6 

11:57:59-12:00:06 

4 

290 

12:02:03-12:03:59 

12 

6 

12:00:50-12:02:55 

5 

200 

12:12:01-12:15:56 

12 

6 

12:12:17-12:16:41 

6 

290 

12:17:48-12:21:29 

12 

6 

12:16:43-12:19:50 

7 

200 

12:27:51-12:31:49 

12 

6 

12:27:00-12:29:35 


TABLE 3.2. Selected Runs of the February 9 Test. 



B-57P 

Aircraft Data 


NOAA Lidar Data 

Run 

No. 

Azimuth 

Angle 

Sampling Time 
(MST) 

Start to End 

PRF 

(Hz) 

Number of 
Pul se 
Averaqe 

Sampling Time 
(MST) 

Stzrt to End 

9 

200 

12:14:06-12:17:45 

12 

48 

12:13:39-12:17:23 

10 

290 

12:19:30-12:23:09 

12 

24 

12:17:45-12:22:13 

11 

200 

12:28:05-12:31:43 

12 

24 

12:28:49-12:30:53 

12 

290 

12:33:25-12:37:09 

12 

24 

12:33:47-12:36:49 
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Figure 3.9. Aircraft-measured wind vector on February 7 flight. 
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Figure 3.11. Evaluation procedure for averaginy B-57B wind for each 
range gate. 
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Figure 3.1 
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2. Comparison of radial mean wind velocity, calculated turbulence 
intensity, and lidar spectral width between aircraft measure- 
ment and lidar measurement on February 7, 1984 (290° azimuth 
assuming horizontal homogeneity). 
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Figure 3.14. Comparison of radial mean wind velocity, calculated turbulence 
intensity, and lidar spectral width between aircraft measure- 
ment and lidar measurement on February 7, 1984 (290° azimuth 
assuming vertical homogeneity). 
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Figure 3.15. Comparison of radial mean wind velocity, calculated turbulence 
intensity, and lidar spectral width between aircraft measure- 
ment and lidar measurement on February 7, 1984 (20° azimuth 
assuming vertical homogeneity). 


40 








corresponding sampling volumes to compare under the assumption of vertical 
homogeneity, when the recorded aircraft coordinates are questionable (due 
to INS drift problem). However, a reasonable guess was made to adjust the 
aircraft trajectories. By comparing Figure 3.12 with Figure 3.14 and 
Figure 3.13 with Figure 3.15, there is no significant difference between 
the two methods of comparison (i.e., assuming either horizontal homogeneity 
or vertical homogeneity). 

Since the mean wind direction on February 7 is 30° from true north, 
the magnitude of wind speed at 20° azimuth is higher than that at 290° 
(comparing Figure 3.12 and 3.13). The negative radial wind speed means 
the wind was away from the lidar site. Differences of radial mean wind 
speed between aircraft measurement and lidar measurement are noticed in 
these figures (Figure 3.12 to 3.15). It was also noticed that the magni- 
tude of aircraft-measured wind was larger than the lidar-measured wind in 
most cases. One possibility is inaccuracy in wind measurements due to 
the Schuler frequency of the aircraft navigation system or dje to varia- 
tion in pulse transmission frequency of the lidar system. Despite the 
difference of mean wind from both measurements, the turbulence intensities 
measured by computing the lidar wind time history for each range gate 
agree very well with aircraft-measured intensity during the interval the 
aircraft is passing through the range gate. The outputs of the second 
moment (lidar spectra width) are also available for the NOAA/WPL lidar 
system. However, the magnitude of spectral width is in the range of 2 to 
3 m/s, which is 4 to 5 times the actual turbulence intensity. It is 
believed that the major contribution of these large spectral width values 
is the broadening of the square pulse itself in the spectrum estimator. 
Careful investigation of the signal process needs to be carried out to 
resolve the correlation between lidar spectral width and turbulence 
intensity. 

Data of February 9 are also examined. Figure 3.16 and 3.17 show the 
comparision of both measurements at different azimuth angle, i.e., 290° 
and 20°, respectively. From both runs of Figure 3.16, the reverse flow 
clearly appeared at slightly below 2000 m. This flow phemonenon only 
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Figure 3.16. Comparison of radial mean wind velocity, calculated 

turbulence intensity, and lidar spectral width between 
aircraft measurement and lidar measurement on February 9, 
1984 (290° azimuth assuming horizontal homogeneity). 
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appears at the direction roughly perpendicular to the mountain range and 
Is not seen in Figure 3.17 which is 20° azimuth or on February 7 when 
the prevailing wind direction Is 30° toward the mountain range. Similar 
flow patterns were detected by B-57B measurements during the turbulent 
flux iasurements portion of the orographic campaign (Theon 1985). 

Turbulence intensities from the February 9 data were also computed 
from aircraft-measured wind and lidar-measured wind. As indicated in 
Table 3.2, the number of pulses averaged are 24 or 48 for Ildar data 
while pulse repetition frequency is 12 Hz. The time interval between two 
sequential outputs is 2 or 4 seconds. A significant loss of turbulence 
information could occur due to long time period averaging. Although the 
spectr-i widths are still higher than computed turbulence Intensities, it 
is interesting to see that the variation of lidar spectral width is 
strongly correlated with the wind shear. As shown in Figure 3.16, a 
strong wind shear occurred slightly above 2000 m altitude for both Runs 
10 and 12 and a larger value of spectral width appeared at the same 
altitude for both runs. 

Turbulence spectra were computed for each of six flight paths (three 
in 290° azimuth and three in 20° azimuth) and at each corresponding range 
gate, assuming horizontal homogeneity. The spectra computed for each 
range gate for the three flights in 290° and three flights in 20° were 
then segment averaged, respectively. Figure 3.18 shows the calculated 
turbulence spectra of some of the range gates for the 290° azimuth paths 
while Figure 3.19 shows results for the 20° azimuth paths. Note that the 
February 7 lidar data are sampled at two times per second resulting in a 
Nyqulst frequency of 1 Hz. The aircraft data, on the other hand, are 
sampled at 40 times per second resulting in a Nyquist frequency of 20 
Hz. The differences of sampling rate and the length of data resulted in 
the spectra falling on different frequency bands; the spectra, nowever, 
do merge at intermediate frequency values forming a relatively continuous 
line and follow the -5/3 slope quite well. Thus, the turbulence measured 
by both instruments displayed highly realistic features of natural turbu- 
lence structure. 
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Figure 3.19. Computed radial turbulence spectra at 20° azimuth path. 
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4.0 CONCLUSIONS 


The results of this study show general agreement between winds 
measured with the Doppler lidar and winds measured with the E-57B instru- 
mented aircraft. In all comparisons the winds along the lidar beam and 
those measured al^ng the aircraft flight path are in agreement within a 
factor of 10. In many cases, the agreement is much better and well 
within ±2 to 4 m/s. It is obvious that exact agreement between wind 
fields cannot be expected for a number of reasons: (1) the aircraft 

never flies exactly the same path as the lidar beam, (2) the lidar 
measurement is averaging the wind along a cylindrical volume element 
'.u0 m in extent whereas the aircraft represents a spatial average based 
on Taylor's hypothesis, (3) the aircraft measurements themselves have 
inaccuracy due to such factors as the Schuler frequency of the aircraft's 
INS, (4) the NOAA lidar was reported at the time of the experiment to have 
variations in pulse transmission frequency which also could induce inaccu- 
racies in wind measurements. In view of all these factors, it is con- 
cluded that the lidar and aircraft measurements essentially show valid 
wind speed measurements. 

Turbulence intensities measured by computing the lidar wind time 
history for each range gate and then calculating the rms value relative 
to the mean agree quite well with the aircraft intensities. The spectral 
width or second moment data from the lidar, however, does not correspond 
well with the aircraft-measured intensities being consistently a factor 
of 2 higher. This difference may be due to the natural pulse itself 
broadening the spectral estimator. 

The computed spectra from lidar measurements fit the spectra com- 
puted from aircraft measurement values very well. Also, they follow an 
approximate -5/3 power law as is expected for turbulence in the atmo- 
sphere. The scatter in the spectral data is quite large but this is to 
be expected because of the small amounts of data used in the statistical 
averaging. The general results of the study, however, suggest that 
turbulence measurements with Doppler lidar systems holds considerable 
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promise. Further work in this area of comparing remote sensed values 
with in situ measurements is needed, however, to fully resolve the 
reliability of 1 idar-measured turbulence values. 
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